Schizophrenia - a parameters' game? 



Anca Radulescu, Applied Mathematics, University of Colorado at Boulder 
UCB 526, Boulder, CO 80309-0526, radulesc@colorado.edu 

February 1, 2008 



Motto: "Ivan concentrated his attention on the cat and saw how the 
strange animal walked up to the platform of an A tram waiting 
at a stop, cheekily pushed off a screaming woman, grasped the 
handrail and offered the conductress a ten-kopeck piece." 



M. Bulgakov - "The Master and Margarita' 



Schizophrenia is a severe, currently incurable, relatively common mental condition. Its symptoms 
are complex and widespread. It structurally and functionally affects cortical and subcortical regions 
involved in cognitive, emotional and motivational aspects of behavior. Its cause is unknown, its 
diagnosis is based on statistical behavior and its treatment is elusive. 

Our paradigm addresses the complexity of schizophrenic symptoms. Building upon recent neural 
vulnerability and limbic dysregulation hypotheses, it offers a mathematical model for the evolution of 
the limbic system under perturbation. Dependence on parameters and the concept of "bifurcation" 
could be the key to understanding the threshold between "normality" and "disease" . 



1 Stress and mental illness 

In the wild, living beings survive by responding to perceived threats with adaptive and appropriate 
changes in their behaviors and physiological states. Besides the species-specific factors, the nature 
of these responses depends on the external environment, but also on the internal physiological and 
emotional conditions. Unfortunately, the neuroendocrine mechanisms that control stress responses 
based on these environments are poorly understood for most animals [72j . in particular for humans. 

Altogether, we are aware of a few physiological mechanisms that allow stressful events to affect 
behavior. The stress response impacts on autonomic functions, such as respiration and heart rate 
variabiltiy [73 . It also alters, through the hypothalamic-pituitary-adrenal axis, the secretion of 
corticotropin-releasing factor, the hormone that coordinates behavior, autonomic and hormonal 
function in response to stress [1] . It subsequently affects the production of multiple other hormones, 
such as adrenocorticotropin hormone and Cortisol. 

It has been generally understood for years that, via its autonomic effects, sustained stress can 
severely affect health contributing to a variety of conditions, among which heart disease, di- 
abetes, growth retardation [105j . decresed immunity [71] and various eating and digestive disor- 
ders Through similarly complex mechanisms, stress is also believed to lead to a number of 
psychiatric disorders, including depression, posttraumatic stress disorder, Alzheimer's disease, and 
other anxiety disorders [55] • Interestingly, whether one is male or female is one of the most impor- 
tant health determinants. This has been partly explained exacly by the underlying sex differences 
in the physiological response to stress, specifically by the fact that estrogen exposure attenuates 



sympathoadrenal responsiveness. It has been hypothesised that these differences have been driven 
in females by strong evolutionary pressure, in the need to protect the fetus from the adverse ef- 
fects of maternal stress responses, in particular excess glucocorticoid exposure [M]. As psychiatric 
stress-associated pathologies will be the object of our present study, we will come back to the idea 
of neural exposure to Cortisol in a later section. 

In relation to behavior and mental illness, stress has been noticed to impair heathy individuals 
and psychiatric patients in sometimes surprisingly similar ways. Consider schizophrenia, for example, 
wich will be the focus of this paper: a severe mental disorder with a heterogeneous set of symptoms 
including paranoia, hallucinations, delusional beliefs, thought disorder, emotional flattening and so- 
cial withdrawal. The illness is relatively common, affecting 1.1% of the population, or over 2.2 million 
people every year in the United States alone, and has a devastating impact on social functioning. 
Studies over the past two decades have established that it is a chronic [9] [52] and neurodegen- 
erative [5] [33] [31] disease, structurally and functionally affecting various cortical and subcortical 
regions involved in cognitive, emotional and motivational aspects of behavior [55] [5] |112| [7]. 

Many cognitive abnormalities [SB] [3S] [TS] associated with this illness, such as delusions and 
hallucinations ^Mj, impaired memory [SU], lowered sensory gating and selective attention [T7] [71] 
are also induced in healthy adults under acute emotional stress [45] |108j [T8] . So the question 
raises itself: Can stress single-handedly cause a mental illness just through its perseverance? And 
if not, then what is the intrinsic detail that makes the dramatic difference between normality and 
pathology? 

Historically, theories have ranged over a large scale, and psychiatry hasn't lacked controversies. 
Some attributed schizophrenic symptoms solely to prenatal stress, or to social and other environ- 
mental factors during the patient's childhood and teenage years |59j . Others made them entirely 
the responsibility of genetic programming [T^] |122j . 

Unfortunately, despite intensive research, our knowledge of the underpinnings of schizophrenia 
is now about as relevant as ever. At present, there is no generically sustainable treatment for most 
serious psychiatric illnesses. Progress has been made in understanding some effects and side-effects 
of medication, but the drugs that are being used may only treat the effects of the disease rather 
than its cause. In fact, the main unanswered questions in current psychiatry concern diagnosis of 
mental illness as much as its treatment. Consensus diagnoses are revised periodically [3^ based upon 
observed behavior rather than actual causes (etiology), which are usually still unknown. Most psy- 
chiatric conditions come with complex, widespread symptoms, accompanied by seemingly unrelated 
cognitive abnormalities and psychosocial consequences. The severity of symptoms cannot be stated 
in reproducible terms, and is therefore left to clinical interpretation, which varies among psychia- 
trists. This increases the potential of misdiagnosis, together with the fact that psychotic or very 
young patients are often neither cooperative nor articulate enough to describe their symptoms. 

This is true in particular for schizophrenia. Recently, an idea which gained ground, under diverse 
forms, has been that the etiology of schizophrenia is based on neural vulnerability and degeneration. 
The theory (first conceptualized in the mid 19th century) attributes mental disorders to a hereditary 
predisposition that reduces the individual psychological threshold towards stimuli [113j , to the point 
where even minor daily stresses will directly trigger psychotic experiences [89] . It has been observed 
that this "vulnerability" (or lack of inhibition in the threat detection mechanism fET ) manifests itself 
as an overt illness only under the impact of stress factors |116J, so that schizophrenic disturbances 
eventually result as an overlap of environmental stress onto the individual's premorbid personality 
component. 

In support to this hypothesis come the well-known relationship between stress and first-break 
psychosis [53j or relapse [11 6J, and the sympathetic upturn that occurs prior to symptom exacerba- 
tion. Prospective data suggest that signs and symptoms (such as elevated autonomic activation [33] 
and clectrodermal activity levels [S3j [35]) prodromal to psychotic first episodes and relapses may 
be present in about 60 percent of patients. Also, it has been observed that stressful life events 
and highly critical attitudes toward the patient in the social environment predict relapse [901. On 
the other hand, research with schizophrenic out-patients has shown that antipsychotic (neurolep- 
tic) medication reduces relapse rates. This protective factor may operate partially by raising the 
threshold in the face of environmental stressors [117j . 

Such first outbreak and relapse predictors are currently being used as clinical indicators for 
schizophrenia, together with more traditional ones, such as paranoia, agitation and sleeplessness [54j . 



2 



The subsequent possibility of pre-symptomatic treatment [79] [80], among other things, motivated a 
more careful investigation of the factors implicated in producing this "vulnerability" , and the attempt 
to determine the conditions that precipitate the plethora of outward clinical manifestations diagnosed 
as "schizophrenia". In this perspective, several vulnerability models have been proposed [13j . which 
we will discuss in some detail in the next section. 

2 The limbic dysregulation hypothesis 

Recent studies have increasingly correlated vulnerability to schizophrenia to volume reductions in 
several limbic areas (amygdala, hippocampus and prefrontal cortex |124j ). In light of the previous 
section, this should not appear surprising, since the limbic system is primarily associated with the 
regulation of emotion and arousal, and is also responsible for integrating the internal and external 
environments via its wide connections with the neocortex [5] , as well as with the autonomic ^ and 
endocrine [7S] systems. 

Dopamine and serotonin abnormalities [63j [123| in schizophrenia constitute today the most es- 
tablished and popular etiological hypothesis (which forms the bases for development of newer antipsy- 
chotics [43] [65]). However, schizophrenia has many neurobiological features suggesting an underlying 
dysregulation of emotional arousal, including limbic [H] |121| , endocrine [100] |104| [120| [115] , and 
autonomic [33] |125| [55] abnormalities. It is possible that the neurotransmitter disfunction may 
be induced by hyperarousal [5S] [S3] [U], making it a consequence of dysregulation, rather than its 
cause. 

Over the years, this variety of abnormalities of schizophrenia (and the relationship among 
themselves and with hyper/hypo-arousal) have spurred research interest. Historically, Gruzelier 
et al. [in] [H] were among the first to hypothesize that limbic abnormalities in the regulation of 
arousal may be an important feature of the disease. Grossberg's model [17] linked schizophrenia to 
pathologies of the amygdala and Hanlon and Sutherland [H] - to prenatal damage to the limbic 
system. Williams et. al. [121, have reported differences between the amplitude of patient and control 
limbic and autonomic responses to different facial expressions. Other studies pointed out the abnor- 
mally high Cortisol levels in schizophrenic patients, [120| |115| . and McEwen |78j has sustainably 
related mental illness to chronic stress and to the corresponding deleterious effects of Cortisol on the 
hippocampus and prefrontal cortex to mental illness. 

In (1989), Nuechterlein and others [50] [SI] [H2] elaborated a "vulnerabihty/stress" hypothe- 
sis of schizophrenia, in which a vulnerability to stress (assumed to be due to cognitive deficits), 
combined with stressful "life-events," leads to first-break or relapse of schizophrenia. Advances in 
understanding the neurobiology of the stress cascade led to a plausible model by which this vul- 
nerability may occur through neurotoxic effects on the hippocampus that may involve synaptic 
remodeling [27] [HS] \TW\ . 

In this context, schizophrenic symptoms may constitute an end-stage of a cyclic and neurode- 
generative process. Recent studies [81] [96] |114j support the theory that the vulnerability to stress 
in schizophrenia is based on a pre-existing hippocampal/prefrontal deficit. Impaired hippocam- 
pal/prefrontal function leads to decreased inhibition of the amygdala, contributing to higher arousal 
levels, even under minor stress. Via the connections of the amygdala with the hypothalamus, the fear 
reaction triggers autonomic and endocrine effects [7D] (such as changes in heart rate variability [53] 
and electrodermal activity [53] ^2j, or increased Cortisol levels [104] ). Excessive Cortisol leads to 
brain neurotoxicity [119j and further hippocampus damage [51] [104j , thus closing the dysregulation 
vicious cycle. The delay in schizophrenia's onset (late teens in males and early 30's in females) is 
consistent with a vicious cycling process, in which the neurodegenerative loop would need sufficient 
time to progress to the point where symptoms become apparent; the already mentioned neuropro- 
tective effect of estrogen in ameliorating the effects of Cortisol on the hippocampus and prefrontal 
cortex [To] [46] would slow but not stop the cycle in females ^67j, accounting for their delayed onset. 

The dynamical analysis in this paper is based on a control system model described by Sotres- 
Bayon et al. p,09], in which limbic regions define a negative feedback loop that regulates arousal. 
The central amygdala forms the main excitatory component of the arousal response |30j . The 
primary inhibitory pathways are the medial prefrontal cortex [II][I1] [SS] [57] [M] |103j and the 
hippocampus [27j |109j . Outputs from the limbic system, via the hypothalamus, provide inputs for 
the endocrine and autonomic nervous systems. In this sense, the model presents a review of the 
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known possible mechanisms for regulation of arousal. In our context, the model explains how limbic 
dysregulation in schizophrenia could lead to its characteristic behavioral features and could also 
cause the endocrine and autonomic abnormalities that so often accompany the illness. 

Section 3 will give a brief review of the known neural pathways that underlie the limbic connection 
assumed in our model. Section 4 will construct and analyze the mathematical model. Section 5 will 
interpret the results of the analysis and will discuss the conclusions in a clinical context. 

3 Connections and pathways 

Over the past decade, significant research has been conducted on the role of the prefrontal cortex, 
the hippocampus and the amygdala in the fear conditioning and extinction. The predominant view 
is that the amygdala is excitatory and the hippocampus and prefrontal cortex are inhibitory [109j . 
More precisely, we believe that the activity of the prefrontal cortex modulates the amygdala fear 
reaction to a stressor. In this section, we will describe in some detail the internal anatomical 
organization and the pathways between the regions involved in the stress-reaction. This will provide 
us with some background and motivation for our mathematical model, although the model itself will 
be much more schematic and will try to avoid detail. 

1. Amygdala. It has been observed, in both human and animal studies, that damage to the 
amygdala prevents the acquisition and expression of fear. It was thereby concluded that the 
amygdala may be the underlying site for fear conditioning and extinction. Amygdala is divided 
into a few physiologically and functionally distinct parts: the lateral amygdala (LA), the cen- 
tral amygdala (CE), the basal nucleus (B) and the intercalated cell mass (ITC). The current 
hypothesized mechanism of the fear reaction, in a very simplified form, is the following: In the 
absence of stimuli, the intra-amygdala connections are suppresing its activation, maintaining 
it at fairly low levels. When an emotionally potent conditioned stimulus is received, it is 
transmitted via thalamic pathways to the LA, then to the CE (either directly or via more 
complex intra-amygdala connections). Finally, the CE has output connections to a set of 
regions that control specific autonomic, endocrine and behavioral responses (autonomic and 
endocrine systems, PFC). The role of B is still controversial. Although there is anatomical [H] 
and physiological [551 evidence that there are strong reciprocal projections of B with the hip- 
pocampus and with the mPFC, B lesions seem to have no effect on fear extinction. It has 
been suggested that the role of B may be to integrate information from the LA, hippocam- 
pus and mPFC, thus being a site of contextual contributions to conditioning. As both the 
hippocampus and the PFC are belived to be crucial in the dynamics of schizophrenia, and as 
contexetual interpretation of threat has been proved to be impaired in schizophrenic patients, 
these interconnections may prove to be of interest to our present study. 

2. Prefrontal cortex. 

Damage to the prefrontal cortex (PFC) is known to generally induce emotional and cognitive 
changes. In fact, these changes seem to be very finely-tuned and region-specific. The PFC 
consists of several functionally distinct subregions, wich include the lateral prefrontal cortex, 
the orbital frontal cortex and the medial prefrontal cortex (mPFC) [SH| [106j [lOlj . The lat- 
eral prefrontal cortex is involved in working memory and executive control functions (such as 
motor control) [82] ■ The orbitofrontal cortex is involved in motivation, reward and emotional 
decision- making [T3]. The mPFC is itself divided into a few subregions: anterior cingulate 
cortex (ACC) and several more ventral areas (infralimbic, prelimbic). The dorsal part of the 
ACC is involved in attention and cognitive control, and the ventral part in emotional regu- 
lation [12]. The functionality of the other subregions hasn't yet been clearly esablished, but 
the predominant view is that neural activity in the mPFC regulates the amygdala-mediated 
fear responses via direct projections to the LA or the ITC, as well as the activity in the 
hippocampus, via projections to CAl (see below). Conversely, experimental studies suggest 
that initiating and sustaining behavior require several types of mPFC modulation, including 
mPFC self-stimulation [SI] [10]. 

3. Hippocampus. The hippocampus is critical in episodic memory cosolidation and for 
aspects of working memory ^Tlj . Unlike the role of the amygdala and PFC in stress processing, 
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which have been confirmed by a wide variety of studies, the potential contribution of the 
hippocampus remains relatively unexplored. 

Structurally, MRI studies [30] found decreased hippocampal volumes in depressed patients 
and correlated the volume loss with the length of the illness. The same volume reduction has 
been observed in schizotypal disorders [37] . This is consistent with the hypothesis that hyper- 
cortisolism could result in hippocampal neurotoxicity in conditions such as bipolar disorder 
and schizophrenia. 

However, although cronic stress has been shown to structurally damage the hippocampus, 
this damage is believed to be restricted to particular subfields [110] [77] , which is possibly not 
sufficient to explain psychotic symptoms. Cerqueira et al. [21j showed that chronic stress may 
also impair working memory and behavioral flexibility indirectly, by affecting not the volume 
or the number of neurons in the hippocampus itself, but rather the synaptic plasticity within 
CAl 66J or of the hippocampus-PFC interactions (see the paragraph below on hippocampus- 
PFC pathways). 

4. Amygdala prefrontal cortex 

Different amygdala nuclei are robustly connected with different regions in the mPFC, sug- 
gesting that the two are functionally coupled. Several studies have shown that the functional 
mPFC activity is inversly related to amygdala activity [S] , and this regulatory interaction is 
believed to be critical for the organism's ability to adapt to change. Although it has been 
proposed that mPFC inhibits activity in the amygdala, the mechanisms of this suppression are 
not yet known. As most mPFC projections to the amygdala are excitatory, it has been pro- 
posed that the inhibition occurs by activation of inhibitory neurons within the amygdala [103j . 
However, based on experimental evidence, a new study [118j suggests a more complex, bidi- 
rectional modulation of fear, in which PL excites amygdala output (via its projections to B) 
and IL inhibits amygdala output (through its projections to LA and ITC). 

It has been argued that dysfunction of the mPFC-amygdala interaction may trigger the emo- 
tional preservation (ususally a hyperactive amygdala and a hypoactive PFC) found in depres- 
sion lOTj, anxiety [29^ and other fear disorders [97] . 

5. Hippocampus prefrontal cortex 

Clinical and experimental studies implicate both hippocampus and PFC in several aspects of 
learning and memory. Not surprisingly, the two units are stongly interconnected and modu- 
late each other's activity in a complex manner. Hippocampal inervation of the PFC is mainly 
excitatory and originates fron the temporal CAl/subiculum region and projects to the pre- 
limbic, medial orbital and infralimbic areas [61j . Conversely, hippocampal memory supression 
is (at least for non-psychiatric populations) under the control of prefrontal regions [5S] . 

Cerqueira et al. [21] explain how stress can influence the integrity of the hippocampus-PFC 
pathway, and thereby explain some of the neurobiological deficits triggered by stress that can't 
be attributed to hippocampal lesions. The study correlated stress exposure with an observed 
volumetric reduction in the upper layers of the mPFC which could not be accounted for by 
neural loss, but rather by dendritic atrophy and retraction of the pyramidal neurons in layers 
H and HI of the mPFC (also see f53^). Although the hippocampus-mPFC pathway was shown 
to be impaired even by a single episode of acute stress |102j , this stress- induced atrophy seems 
to be reversible [55] . 

6. Amygdala — hippocampus The amygdala impact on the hyppocampus is best represented 
not by neural pathways, but by the indirect autonomic and endocrine effects initiated in the 
amygdala in response to stress, which lead to hippocampus impairment and functional reduc- 
tion (as described above) . Conversly however, studies such as [26] have opened the possibility 
that hippocampal projections to the B might be important for contextual contributions in 
fear extinction. 
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Figure 1: Two schematic models of the amygdala-hippocampus-mPFC interactions, as illus- 
trated by Sotres-Bayon et al. |109j . 



4 The mathematical model 

Our theoretical, simplified model quantifies the amygdala-mPFC mutual regulation in a way which 
can be studied and understood analytically. The pathways between the two regions, as well as the 
self-modulation within each region are represented by linear terms, while the indirect hippocampus- 
modulated influences are expressed as nonlinearities. As we don't consider bilaterality, our model is 
constructed as a 2-dimensional dynamical system in which the variables a = a{t) and p = p{t) are 
levels of activation of amygdala and prefrontal cortex, respectively. (In a data-driven model, one 
may think of these variables as BOLD percent signal change.) 

The strengths of the interactions are tuned differently for different individuals. In our model, this 
tuning is quantified by a set of parameters, so that the time-evolution of the system, and ultimately 
its asymptotic behavior, depend on the choice of these parameter values. While making no claim 
to illustrate exactly the complex fear reaction, the model should rather be seen as a metafor of 
the brain undergoing stress, supporting the limbic dysregulation hypothesis. This has no imediate 
practical value. However, if explored further, the idea may turn out to be clinically priceless, as it 
suggests ways in which a more quantitative approach would be helpful to the field of psychiatry. 
New clinical paradigms, could be developped to test and use this hypothesis, as we will later explain. 

In the next sections, we will show how for different values of the parameters, the system can ex- 
hibit dramatically different time evolutions: some corresponding to normal physiology and behavior, 
some corresponding to typical symptoms of schizophrenia. 

4.1 A linear model. Pros and cons. 

Let's start with the strong (and, as we will see, unreasonable) assumption that the amygdala-mPFC 
dynamic is linear, and criticize the flows inherent to this oversimplification. One would not expect, 
in general, for complex phenomena in the brain to behave in a linear fashion. It is a well-known 
mathematical fact that linear systems don't exhibit any "interesting" behavior. As a rule of thumb 
an isolated fixed point is either a global attractor or repeller, a center or a saddle point (depending 
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on the nature of the Jacobian's eigenvalues at that point). Moreover, hnear systems have no hmit 
cycles, due to the Poincare-Bendixon theorem. The only circumstance in which cycles exist is if the 
system has a "center" , i.e. a fixed point surrounded by accumulating cycles, neither attracting nor 
repelling them. (Mathematically, centers happen when the Jacobian has two conjugate imaginary 
roots.) 

In this context, if the parameters of the system change, the behavior of a fixed point may change 
from being a global attractor to being a repeller, going through a transitional stage (bifurcation). 
The bifurcation corresponds to a critical set of the parameters for which the fixed point is a center, 
and thus represents a point of sudden transition in the dynamics (see Figure 1). As one would want 
to believe that brain processes have somewhat more fail-safe mechanisms, this is a feature of the 
linear model that we will try to correct later. 

As an illustration, we construct a very simple and plausible linear model for the basic interaction 
between the amygdala and mPFC. This toy version of the model will leave out the autonomic and 
endocrine feedback loops, and will therefore neglect their influence on the hippocampus function. 
All other region interactions will be assumed linear. 

d = —/lift — kip + / — 7i-ff 
p = k2a + iJ,2P + J2H 

where /, i/, /i, fci, A:2, 71, 72 > 0. 

The amygdala activation a is driven by four terms: the input / > (corresponding to the 
background environmental stimuli), the self- inhibition — /iia (the amygdala "resilience to stress") 
and the prefrontal and hippocampal modulations ~kip and 71 (the hippocampus is assumed to 
have constant activation H, since it receives no feedback). In the absence of the PFC term, the 

asymptotic equilibrium for a would be a = . The PFC activation p is driven by: the 

^1 . 

amygdala excitatory input k2a, a hippocampal excitatory modulation J2H and a self-excitation 

t^2P- 




-3-2-1012 -3-2-10123 -3-2-10123 



Figure 2: Phase portraits of the linear model for fi2 = 1, A;i = 2, A;2 = 4 and different values 
of the parameter The fixed point (a* ,p*) is a global attractor for fii = 1.1 (left), a center 
for /ii = 1 (middle) and a repeller for fii = 0.9. 

In a "brain" with a well-balanced feedback, this should euristically work as follows: If the 
amygdala activation level is high, its excitatory effect on the mPFC will drive p up, which in turn 
will inhibit a and make the system converge towards an equilibrium. In a "brain" with an overactive 
amygdala (low /ii value), convergence fails to happen. 

The Jacobian matrix: 




has determinant A = kik2 — A*iAi2- We will assume that A = kik2 — /ii/X2 > 0, i.e. we work in a 
regime where the interconnections prevail over self-modulations in the two regions. 
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Then the system has always a fixed point (equihbrium) at 



ia,p)- A ' A ) 

The equilibrium depends on the outside stress level and on the hippocampus activation, as well 
as on the mutual and self excitation and inhibition parameters. The fixed point is globally attracting 
if fii > fi2, i-e. if the amygdala self- inhibition is strong enough to exceed the PFC self-excitation. 
However, as /ii decreases, the stability of the fixed point changes, and for fii < fi2 it becomes a 
global repeller. A system operating under such values of the parameters (/ii < 1^2) will have all 
trajectories pushed away from the repeller and never stabilize. It is hard to believe that the brain 
operates with switch between stable and unstable behavior so suddenly. 

So where could the nonlinearity come from and what does it do? 



4.2 The nonlinear model 

Nonlinearity in a two-dimensional system introduces an interesting characteristic feature: possible 
existence of limit-cycles. Consider the general form of a 2-dimensional nonlinear system with a fixed 
point at the origin: 

\y J \y ) ^'^^ V d2ix + d22y + N2ix,y) J 

where D is the Jacobian matrix and N is the nonlinear part, both parameters-dependent. 

The eigenvalues Ai and A2 of D may take real or complex values that satisfy A1-I-A2 = trace(D) = 
T and A1A2 = A = det(Z3). In general, if A > 0, than the origin is either an attracting node or 
spiral (in case t < 0) or a repelling node or spiral (in case r > 0). At the parameter values where 
T = Ai -|- A2 = 0, the system exhibits a bifurcation. In particular, if Ai and A2 — Ai are imaginary, 
then we may have a Hopf bifurcation. The way the local dynamics of the system changes at a Hopf 
bifurcation is described by the Lyapunov number, which depends on both the linear part D and the 
nonlinear part N, as described below t93|. If A'^i,iV2: R. are analytical, with expansions: 

Ni{x, y) = ^ aijx'y^ = {a2ox^ -f auxy + ao2y^) + (aaoa;^ + a2ix'^y + b^xy^ + boay^) + ... 

i+3>2 

N2{x, y) = ^ bijx'y^ = {b2ox'^ + bnxy + 5o2y^) + (6302:^ + 6212^^^ + bi2xy'^ + 6o3Z/^) + •■• 

i+j>2 

then the Lyapunov number: 



= TT, — rT7?{'^ii^2i(aii + '^11^02 + ao2&ii) + -Diidi2(6ii -I- a2o&ii -t- aii6o2) 
2ai2A-'/^ 

+ c'2i(«ii«02 + 2002602) - 2(iii(i2i(6o2 " £120002) - 2c?iidi2(a2o - 620^02) 

- ^12(2020620 + 611620) -I- (^12^21 - 2dii)(6ii6o2 - aiia2o) 

- (dii + rfi2rf2i)[3(d2i6o3 - ^12030) + 2(iii(a2i + 612) -I- (^21012 - ^12621)]} 

If cr 7^ 0, a Hopf bifurcation occurs at the critical value of the parameters where r = 0. More 
precisely: 

(1) If tJ < 0, the origin is attracting for r < 0, and the system has no limit cycle. For r > 0, the 
origin becomes repelling, but a circular stable limit cycle forms around it, whose radius increases 
with T. 

(2) If (T > 0, the system has a unique unstable limit cycle that surrounds the stable origin for 
T < 0. The radius of the cycle decreases with t. At t = 0, the unstable cycle collapses into the 
origin, making it unstable for t > 0. 
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TrlD><0 



Tr(D)>0 



Figure 3: When a < 0, the system exhibits a supercritical Hopf bifurcation. 




a: 



Figure 4: When a > 0, the system exhibits a subcritical Hopf bifurcation. 



The change in dynamics is still very sudden at a Hopf bifurcation (even though the parameters 
change smoothly), but has much more subtle implications. To interpret this phenomenon in our 
clinical context, let's fist construct as before a more particular example of a nonlinear model, as an 
extension of our linear toy-model. 

Staying faithful to our original assumption, we include the nonlinearities as the contribution of 
the autonomic and endocrine response to stress. As mentioned in the introduction, experimental 
studies showed that high levels of Cortisol due to the amygdala stress-reaction have a detrimental 
effect on the PFC, in two known distinct ways: its neurotoxic effects on hippocampal cells and 
its suppression of synaptic function between hippocampus and the receptive PFC areas. The new 
system, including these nonlinear contributions, will take the form: 



Here the term —/(a) signifies the structurally detrimental effect that amygdala overactivity and 
the subsequent hypercortisolemia has directly on the hippocampus. The term g(a,p) refers to the 
synaptic remodeling induced by high Cortisol on the PFC afferents from the hippocampus, and hence 
it only affects p. 

Under the change of variables x = a — a* and y = p — p*, the system becomes: 

X = nix - kiy + 7i/i(a;) 

y = k2X- ii2y - l2{h{x) - j{x, y)) 

where h{x) = f{a) = f{x + a*) and j(a;,y) = g{a,p) = g{x + a* ,y + p*). 

To fix our ideas, we will take both terms to be of the simplest possible nonlinear forms: h{x) = jx'^ 
to reflect the decease in H activation and j{x, y) = dxy, for the suppression of Hebbian-like synapses 



a = -fiia - kip + I - 71 (7? - /(a)) 
p = k2a- fi2P + 12{H - f{a) + g{a,p)) 
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in mPFC. Recall that, in our original system, this corresponds to /(a) = 7(0 — a*)'^ and g{a,p) = 
Sia-a*){p-p*). 

The system has a fixed point at {x*,y*) = (0,0), and exhibits a Hopf bifurcation at fii = 
with Lyapunov number 

Stt 

^ = 2A3/2 ^^^"^ + ^'^'^^^ [M2('572 - 771) + ^1727] 

(5 k 

Applying the theoretical results we discussed above, we conclude that: if — ! < — (i.e., 

7 Ml 72 

(T < 0), then the system has a supercritical Hopf bifurcation at IJ.1 — H2, and if — I — 3- > — (i.e., 

7 A^i 72 

cr > 0) then the system has a subcritical Hopf bifurcation at fJ-i — IJ.2- 



5 Discussion 

This theoretical model is more a "phylosophical" than "physiological" illustration of the limbic 
dysregulation and neurotoxicity hypotheses. It presents the working brain in a light that permits 
interpretation of its "stress vulnerability" and "hippocampus deficit" as parameters that vary contin- 
uously, determining its regulation and function. The focus of this interpretation is on the idea that, 
although these parameters change smoothly over a whole continuum of possible values, there are 
critical/threshold values, which, when passed, could suddenly and completely change the system's 
dynamics. 

The two parameters on which we focus, whose tuning determines the behavior of the system, are 
the amygdala sensitivity to stress, represented by /ii and the more diffuse hippocampal/prefrontal 
vulnerability to Cortisol neurotoxicity, represented by the Lyapunov number cr. Indeed, staying 
faithful to the stress/ vulnerability hypothesis of schizophrenia, our interpretation will regard a as 
the "disease-quantifying" parameter: negative values of a correspond to normal limbic regulation, 
while positive values of a quantify risk for developing schizophrenia, and, in more advanced stages, 
severity of the illness. On the other hand, larger values of /ii correspond to a more stress-resilient 
amygdala, while smaller values of /ii signify a more stress-reactive amygdala. We could think of 
this parameter as quantifying the amygdala responsiveness to stress, which in literature has been 
related to mental conditions such as depression, or anxiety disorders [109j . but which is not the 
signature of schizophrenia. We can verify whether this paradigm is clinically plausible by testing 
what happens if we apply a brief stress increase to the system (which in real life may come in the 
form of a taumatic event). We quantify the burst of stress by boosting the amygdala to a high initial 
state. We observe wheter the system returns to homeostasis by checking if the respective trajectory 
eventually stabilizes. 

Suppose cr < 0. For high amygdala resilience /ii > fi2, all trajectories converge to the global 
attractor, so the initial condition is irrelevant: the time evolutions stabilize after any stimulus, if 
sufficient time is allowed to pass. For /zi < ^2, the situation is changed by the formation of an 
attracting limit cycle(Figure 4a). After a short stress burst, the duo amygdala-PFC slowly stabilizes 
towards the cycle. Note that, although dampened in time, the system continues to oscillate in both 
cases. We will return to this idea later. The memorable feature of the cr < regime is that, although 
stability of the fixed point changes at the bifurcation fii = ^2, the role of the attractor is assumed 
by a limit cycle. The fact that the amygdala-PFC pair exhibits in some people wider oscillations 
that don't seem to dampen in time could be a mark of low amydgala self-inhibition. 

When cr > 0, the situation changes completely. In the regime Hi > ^2, the system has a locally 
attracting fixed point, surrounded by a repelling cycle, whose radius gets smaller with smaller ^1 
and with larger cr. The basin of attraction of the stable point is the interior of the cycle: any 
initial state inside this basin will converge in time towards the fixed point, and any initial state 
outside will spiral out to infinity. When /zi = /i2, the cycle disappears, so when /ii < ^2 the 
fixed point is globally repelling (Figure 5b). The behavior of the model in the positive a regime 
is representative for schizophrenic dysregulation. A stimulus may elevate amygdala to a value 
which places the corresponding state outside of the attraction basin, preventing convergence from 
ever happening. If alowed to follow its natural evolution, the trajectory would perform larger and 
larger oscillations, corresponding to the cyclic psychotic behavior observed in patients. At some 
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Figure 5: We illustrate the dynamics of the nonlinear model for a particular example of 
parameter values. Fixed parameters: fi2 = = 2,^2 = 4, 71 = 2,72 = 0.5,7 = 1- 

A. The system exhibits a supercritical Hopf bifurcation for 5 = 1. The phase portraits 
show the local dynamics around the origin for jii = 0.9 (panel A.l) and = 1.1 (panel 
A. 2). As III crosses the critical value ^2 = 1; the attracting origin becomes repelling and 
surrounded by an attracting limit-cycle. B. The system exhibits a subcritical Hopf bifurcation 
for = 3. As crosses the critical value /xi = 1, the repelling cycle surrounding the 
attracting origin (fii = 0.9, shown in B.l) collapses into the origin and changes it to a 
global repeller (fii = 1.1, shown in B.2). 



point during this evolution, the patient may enter clinical treatment; antipsychotic medication may 
succeed to temporarily alter this time-course. It is interesting to note that, as we increase the value 
of (J, the attraction basin shrinks around the attracting point, so that even weak stressful stimuli may 
push the state outside the attraction basin. This makes the Lyapunov number a a good quantifier 
of the risk and severity of the disease: the larger tr, the more likely it is for the time-evolution to be 
thrown outside of the convergence range even by small perturbations. This relates to the fact that 
highly vulnerable individuals may develop psychotic behavior even after common daily stress that 
may appear benign to others. 

We want to conclude the analysis of the phenomena described above with a few time-frame 
considerations. Our model works in parallel with behavioral research that relates stress with first 
outbreak [S3] and relapse [116] . The time- frame that the model addresses is thereby of the order of 
days, even weeks, allowing possible psychotic behavior to develop after a stressful event. Cycles and 
oscillations in brain activity, as well as regulation and return to homeostasis have been already the 
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subject of a wide variety of studies, addressing very different time scales. Elecrophysiology studies 
revealed high frequency oscillations in neuronal brain activity, and correlated synchronization of 
theta rhythms in the amygdalo-hippocampal pathways with retrieval of conditioned fear [25j [92j . 
It is also known that that the brain has a circadian rhythm, so that its activity oscillates according 
to a daily pattern [50j . Additionally, recent imaging studies focus on the dynamics of certain 
regions of interest in response to visual stimuli, such as facial expressions. (See Figure 6, which 
illustrates comparatively the behavior of the amygdala and a prefrontal area, in control subjects and 
schizophrenic patients.) Although our model does not address such short-scale phenomena, these 
findings better describe its place as just one level in a general complex picture. 
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Figure 6: A. fMRI data on amygdala and BA^S activation under visual stimulation with 
Pictures of Facial Affect. The data has been averaged over all stimulation blocks, as well as 
over the control and population samples, respectively (samples sizes N = 11) [99j . 



Further testing of the model should be willing to compromise precision for duration. Firstly, the 
design does not need to focus on identifying the exact brain areas involved in the process, but rather 
in estimating the average activation in generic areas such as amygdala, hippocampus and mPFC. 
Secondly, these activations need to be measured over an extended time-period, allowing the feed-back 
cycle to repeatedly loop around, so that the long-term dynamics could become apparent. In this 
context, EEG readings would be ideal, as they can be taken at length, without drastically restricting 
the subject's normal routine. To create the required set-up and take the measures following a stress 
burst, the testing should be ideally planned for immediately after a scheduled stressful event such as 
an exam, or first time ski-diving (see http://www.bme. sunysb.edu/people/lstrey/research.htm^2). 
Another variable of interest that can be very conveniently measured during a long-term routine 
is salivary Cortisol, which is believed to reflect stress-related hypecortisolcmia. In our model this 
variable is one of the factors determining the system's nonlinear behavior, so observing its evolution 
would help us understand these nonlinearities better. To the best of my knowledge, these testing 
paradigms are - if not already being used in other tests - then fairly easy to design and reproduce. 



6 Conclusions 

A natural question to addresses is the significance of this model. Theoretically, the model shows 
how two different clinical systems (with very similar underlying rules, and only slightly different 
parameter values) can exhibit drastically different long term behavior if started under the same 
initial conditions. The literature talks about the "continuum" of human behavior and the practical 
difficulties of establishing a normality/pathology threshold. Such a bifurcation could constitute the 
needed threshold for clinical evaluations. Practically, if the model proves to be valid, both diagnosis 
of illness and quantification of its severity can be achieved by calculating the Lyapunov number of a 
system constructed from clinical measures. As computer capabilities have sky-rocketed over the past 
few decades, fast and accurate algorithms [33] [23] have become available to compute such system 
invariants from time-series. 

On a slightly different note, our model supports the idea that the dynamics of a diseased system is 
not driven randomly, but is rather only apparently random due to its complicated behavior over short 
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time periods. This idea is very important for clinical treatment, as it suggests that the deterministic 
behavior of a system can be changed by proper tuning of the parameters. In general, the drugs 
currently in use only temporarily alter the time-course of the illness, which typically relapses when 
interrupting medication. It is possible that medication options could be improved by exploring how 
drugs can change the parameter values to permanently alter the system and its long-term behavior. 

A more sophisticated, higher-dimensional model, although preferable and possible, is a different 
and much harder task to undertake. As discussed before, the exact physiological underpinnings of 
the brain interactions are not yet fully known from experimental studies, making an educated choice 
of parameters very difficult. Furthermore, even the theoretical study of the evolution of a system's 
dynamics under perturbations is an ongoing mathematical problem. However, we expect that, as 
the theoretical and experimental knowledge progress, so will our ability to understand and model 
such complex phenomena. 



Epilogue 

While walking in a dark swampy place one night, I had the distinct feeling that something was 
following me, with insane reddish eyes, ready to leap. The darkness, the rotten smell, the chill in the 
air, all made my skin wrinkle in little goose-bumps. A stick cracked under my foot, and my heart 
skipped a beat. But I did not run off, neither did my heart stop. A split second later I skipped 
back to the reality of my short stroll in the vicinity of a very populated area, which never harbored 
red-eyed beasts. I giggled at the memory of all the horror Hollywood movies that lived involuntarily 
on that brief adrenaline peak. The moment was gone and the heart was beating all its strokes again. 
But I will never know what just happened in that brief moment of panic. How many more years of 
research arc needed to shed light onto something as seemingly simple, and when will we be able to 
understand why such a casual moment sometimes conquers people and takes over their lives, keeping 
them is the constant terror of some red-eyed beast. 
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